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Abstract. The long time behaviour of the survival probability of initial state and its 
dependence on the initial states are considered, for the one dimensional free quantum 
particle. We derive the asymptotic expansion of the time evolution operator at long 
times, in terms of the integral operators. This enables us to obtain the asymptotic 
formula for the survival probability of the initial state ip{x), which is assumed to 
decrease sufficiently rapidly at large \x\. We then show that the behaviour of the 
survival probability at long times is determined by that of the initial state ip at zero 
momentum k = 0. Indeed, it is proved that the survival probability can exhibit the 
various power-decays like t~ 2m ~ l for an arbitrary non-negative integers m as t — > oo, 
corresponding to the initial states with the condition tp(k) = 0{k m ) as k — > 0. 
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1. Introduction 

The decaying quantum systems such as an a-decaying nucleus are often described by 
the survival probability of initial state, which is the probability to find the initial state 
in the state at a later time. One of the remarkable properties of the survival probability 
is its power decay law at long times. This is a mathematically predicted nature for the 
systems which possess the continuous energy spectrum bounded from below [|lj. It is 
actually seen in many models (see, e.g. || ^ |4[] and the references therein). On the 
other hand, there still remains the difficulty in observing such power decays in the real 
experiments ||, ^, |7|. Hence, the further theoretical and experimental investigations of 
the power decay law are much required. 

One of the fundamental and important models that exhibits the power decay law 
is the free-particle system, from which we can gain many insights into the dynamics 
of the survival probability and also the spatial wave packet. Recently, another aspect 
in the power-law regime at long times is revealed for the one dimensional free particle 
system, in connection with the initial states. As we known, the spatial wave packets 
for this system are expected to decay like t" 1 ^ 2 at long times t. We are assured of 
such a decay for the Gaussian wave packet. However, it is not necessarily valid for 
an arbitrary initial state. In fact, if we take a spatial power-law wave packet as an 
initial state, the "anomalous decay" of its maximum can occur with the form t~ a / 2 
(1/2 < a < 1) [5|, r0|. This is obviously slower than the well-known t -1 / 2 decay. 



The faster decay than t 1//2 is also studied for the initial wave packets which vanish at 
zero-momentum, in association with the dwell time |TI| and the time operator (see [|T^ 



and Appendix A). The latter refers to the time evolution of the survival probability. 
Therefore, the asymptotic decay form of the wave packet for the free particle system 
depends on the initial states in a considerable way. However, we seem still not to get 
a clear perspective of this new aspect of the power decay law. Our aim in the present 
paper is to find the strict condition of the initial states for the various power decays in 
the one dimensional free particle system, and to clarify the underlying mechanism for 
such decays. In particular, we restrict ourselves to the survival probability which will 
show the various power decays as same as the spatial wave packet. 

To this end, we introduce a systematic approach which consists of the two 
procedures. We first derive in section |3| the asymptotic expansion of the time evolution 
operator exp(— UHq) as t — > oo, where Hq is the free Hamiltonian for the one dimensional 
free particle system. It is formally written as 

oo 

exp(-ittfo) = tT 1 £(-ir ^(j + 1/2)(U)-^ 2 G 2j , (1.1) 

where T(z) is the gamma function, and G2/S are some integral operators. This kind 
of the asymptotic expansion was already developed by Rauch [113], Jensen and Kato 
fHf , and Murata with the detailed analyses (see also a recent comment by Amrein 



|i6|). Their results concern the time evolution operator for the systems with a short- 



range potential V(x). The asymptotic expansion (|1 . 1|) for the free particle system is 
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not evaluated by the authors mentioned above, however, it can be achieved without 
difficulty, following their technique for the potential system. The survival probability of 
the initial state if) is defined by the square modulus of the survival amplitude of if), that 
is 

(if>,exp(-itH )if)). 

We see from Q) that (if), exp(-itH )ip) = 0(iT 1/2 ) only if (if),G ip) ^ 0. In other 
words, if the following condition, 

(if),G 2j iP) = 0, j = 0,l,...,m-l, (1.2) 

holds for some integer m, we can obtain another asymptotic decay form of 
(ip,exp(-itH )ip), t"" 1 " 1 / 2 , faster than t' 1 ' 2 . Thus, our remaining procedure is to 



interpret the condition ( 1.2j) as the behaviour of the initial state ip in, e.g., the position or 
momentum space. This is achieved in section f|, and we will finally obtain a remarkable 
conclusion: if the initial state if) behaves like if)(k) = 0(k m ) at zero momentum k = 0, 
then 

|^,exp(-i^ )^)| 2 = Tr-H-^-^Cm+l/^^^G^r+oCr 2 *"- 1 )^!^) 

as t — > oo, where if) is assumed to decrease sufficiently rapidly at large \x\, but otherwise 
arbitrary. Hence, the behaviour of the initial ip at zero momentum, i.e., ip(k) = 0(k m ), 
completely characterizes the asymptotic decay form of the survival probability of if). 
This fact is also expected from the results on the decays faster than t _1//2 for the one 
dimensional free particle system JTTj, |T2 |. 



The organization of the paper is as follows. We first consider in section |2| the 
asymptotic behaviour of the free resolvent (H — z)~ x at small and large energies. This 
study is necessary to the proof of Theorem |3j] in section ||, where the asymptotic 



expansion (|1 . 1| ) is derived. The derivation essentially follows the method used by Jensen 
and Kato | l4fl . To derive the asymptotic formula (|1.3|) for the survival probability, it is 
enough to derive that for the survival amplitude. The latter is accomplished in Theorem 
in section [|. Concluding remarks are given in section |5]. 



2. The free resolvent in one dimension 

We here define the free Hamiltonian in one dimension Hq := P 2 , where P is the 
momentum operator defined by P := —\D X , D x being the differential operator on 
L 2 (R). Then, the free resolvent Rq(z) := (H — z) -1 is explicitly represented as an 
integral operator on £ 2 (R) |T7 



(Ra(z)if))(x) = --A-r I exp(i// 2 |x - y\)if)(y)dy, (2.1) 



R 



21Z 1 / 2 

for all ip E L 2 (R), where z belongs to C + := {z E C | Im z > 0}, and Im z 1 ^ 2 > 0. The 
resolvent Rq(z) is analytic in z. We, however, intend to regard it as an operator which 
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maps L 2 ' S (R) to L 2,_S '(R) for positive s and s'. Here L 2,S (R), defined for an arbitrary 
real s, is the weighted L 2 -space with the norm 

-j 1/2 

i + x 2 ywx)\ 2 &x 

R 

where ip G L 2,S (R) means \\ip\\ s < oo. For positive s and s', the relation L 2,S (R) C 
L 2 (R) C L 2,_S '(R) holds. In addition, we denote by B(s, — s') the Banach space of the 
bounded operators M from L 2,S (R) to L 2, ~ S '(R), with the norm 

||M|| a ,_y := sup 

V>ei 2 - S (R), V^o 

Notice that M G B(s, — s') means the finiteness of its norm ||M|| S a / < oo. The reason 
for this kind of preparation for the free resolvent will be clear in the last part of the 
next lemma. 

Lemma 2.1 : For s,s' > 1/2, Rq{z) belongs to B(s, — s'), and is continuously 
extended to C + \ {0} where C + the closure of C + . 



Proof : We have an estimation 

\\Ro(zW_ s ,= [ (l + x 2 y s ' 

JR 

2 
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2 



. r 1/2 j^exp(iz 1/2 \x - y\)ip{y)dy dx 

for all ip G L 2,S (R) and z G C+. This result clearly holds for z G C + \ {0} and the last 
part of the statement is also proved. □ 

We use the same symbol Rq(z) for the extension of Rq(z) to C + \ {0}. The free 
resolvent Rq{z) is formally expanded around z = 0, 

oo 

Ro(z) = EO^^i, (2-2) 

j=0 

where Gj (j = 0, 1, . . .) is an integral operator acting on the suitable vectors ■?/>, 

(G^)(x) := ~ [ \x- yH{y)dy. (2.3) 



2j\ Jr 

Lemma 2.2 : The integral operator Gj is a Hilbert- Schmidt operator that belongs to 
B(s, -s') with s, s' > j + 1/2. 

Proof : Note that the statement in the lemma is equivalent to (1 + x 2 )~ s '/ 2 Gj(l + 
y2^-s/2 being a Hilbert-Schmidt operator on L 2 (R). The latter is easily seen from the 
relation 

£ y\ 2 ^ f f I x 1 2 ^ ~~\~ I y I 2 ^ 

dxdy <2 2j / — — — - — ^-dxdy < oo 



/rJr (1 + x 2 ) s '(i + 2/ 2 ) s * ~ 7r7r (1 + x 2 ) s '(l + y 2 ) s 
for every s, s' > j + 1/2. □ 

The validity of the formal expansion (2.2) is ensured at small energies, in the 
following sense. 
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Lemma 2.3 : Let k = 0, 1, . . .. If Rq(z) is approximated by a finite series in ( |2.2| ) 
up to j = k, the remainder is o{\z\^ k ~ ls> ' 2 ) as \z\ — > 0, in the norm of B(s, — s') with 
s, s' > k+ 1/2. In the same sense, (|2.2|) can be differentiated in z G C + \{0} any number 
of times for appropriate s and s' , that is, the r-th derivative in z of the approximating 
finite series is equal to (d r /dz r )R (z) U P t° an error of o{\z\( k ~ 1 ^ 2 ~ r ) in the norm of 
B(s, -s') with s,s' > k + r + 1/2. 

Proof : We first consider the case of k — 0. Suppose that r is a non-negative 
integer, s, s' > r + 1/2, and ip G L 2,S (R). Then it follows that 



dz r 



dz r 



[1 + x 



2\-s' 



R 



d r 

r dz r 



< 



X 



A\z 



R JR 

l-2r 



^- I ^(exp(iz / Ix-j/l) 
(l+x 2 )-'(l + 2/ 2 )- s 

1/2 I 



dx 



exp{iz 



x 



|z| m |x - y| 



dzdy(2.4) 



where A and A m ^ m i are positive constants, and m and m' are integers, satisfying 
— 1 — 2r < m < —1 — r and < m' < 2r, respectively. When r = 0, there is no 
contribution from the summation J2mm'- By the dominated convergence theorem, we 
see that (p.4[) divided by | ^ | 1 2r ' goes to as z — » 0. In the same way, for fc = 1, 2, . . ., 
we have 

d r R (z) 



dz r 



i=o <l - 



1 + z 



2\-s' 



R 



X 



d^ 

r dz r 



(iz 1/2 ) fc 1 |a: — 2/| 
2(fc- 1)! 



t 1 " 1 exp(iz 1/2 \x - ?/|(l -*))di 



ip(y)<iy 



dx 



< 



\2(k-l) 



[1 + x 2 



R JR 



2\-s 



fc-l-2r 



a? -y 



2k 



X 



t fc 1 exp(iz 



1/2 



x — y|(l -t))dt- - 



|z| m |a; - j/| 



dxdy 



(2-5) 



for s, s' > + r + 1/2, where -B and £? m>m / are positive constants, and m and m' are 
integers, satisfying k — 1 — 2r < m < k — 1 — r and 2k < m' < 2(k + r) , respectively. 
Taking the limit z — > 0, one can see that (|2.5| ) divided by |z| fc_1_2r goes to 0. This 
completes the proof of the lemma. □ 

On the other hand, we also have the following lemma with respect to the asymptotic 
behaviour of Rq(z) at large energies. 
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Lemma 2.4 : Let k = 0,1,... and s,s' > k + 1/2. Then Ro(z) is k-times 
differentiable in z G C + \ {0}, in B(s, —s'), and it behaves like 
d r R (z) 



dz r 



0(\z 



-(r+l)/2> 



0,1,. ..,k, 



as \z\ -too in the norm o/B(s, — s') 



Proof : Suppose that s, s' > + 1/2 and ip G L ' S (R), we have 

2 2 



< 



4 



R dx (i + x 2 y 



R 1 



y 



2\s 



X 



D| 2 |- 1 -'-|x- 3 /| ar + £ D m>m ,\z\ m \x-y\ 



< oo, 



where D and D m m > are positive constants, and m and m' are integers, satisfying 
— 1 — 2r < m < — 1 — r and < m' < 2r, respectively. Then, the right-hand side 



is 0(\z\ 



as \z\ 



oo. □ 



3. Asymptotic expansion of the time evolution operator 

In order to derive the asymptotic expansion of exp(— UHq) in (|1 . 1| ) , we first define the 
spectral density denoted by E'(X) := (27ri) _1 (_R (A) — -Ro(A)) for all A > 0, where 

(R (\m(x) :=1 J^J^dy exp(-iA 1 / 2 |x - y\)^{y), 

for every if> G L 2 ' S (R) with s > 1/2. The operator Rq(\) is considered as the limit 
of i?o(A + ie) in e | 0. E'(X) clearly belongs to B(s, — s') with s, s' > 1/2, and it has 



the same properties as Rq(X) described in Lemmas pTT| , [2~^ , and |2l4] . Substituting the 
expansion ( |2.2|) and the corresponding one of Rq(X) into E'(X), we have 

n 

E'(A) = 7T- 1 ^2(-iy-^X j -^ 2 G 2j + F n (X), (3.1) 
i=o 

where F n (X) is the remainder. It should be noted that there are no integer powers in A. 
We next focus our attention on the following formula 

exp(-f£# ) = lim lim [* E\X)e~ itx dX, (3.2) 

R^oo rj.0 J r 

valid in B(s, — s') for s,s' > 1/2. The integration in the above can be also regarded 
as the complex integral of (2iri)^ 1 R (z)e~ ltz with the contour enclosing the spectrum of 
H , i.e., [0, oo). This formula is shown in Appendix B. Then, the asymptotic expansion 
of exp(— itHo) at large t is obtained from the formula (|3.2|) together with the expansion 
(|3.1|). To be precise, we can show 



Theorem 3.1 : Let n = 0, 1, . . . and s, s' > max{3n + 3/2, 5/2}. Then it follows 
that 

n 

exp(-ittfo) = vr- 1 EC-^rO' + l/2)(itp- 1/2 G 2j + o^- 1 / 2 ) (3.3) 

j=0 

as t — > oo, in the norm o/B(s, —s'). 
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-itH ) at large times is only determined 
at small energies. Our proof follows the 
and it is given in Appendix C. 

by using the "generalized" zero-energy 
Since Go = — ^(vo,- )v 9 o from (|2.3|), we 



)<po + o(r 1/2 ) 



t 



oo. 



(3.4) 



Note that the asymptotic form of exp( 
by the behaviour of the free resolvent Ro{z 
procedure proposed by Jensen and Kato [|14| 

It is interesting to rewrite the formula 
eigenfunction of H , i.e., (fo(x) := (27r)~ 1//2 
have an alternative expression of ( |3.3| ) as 

exp(-iti7 ) = ^(it)- 1 ' 2 (if. 

Note that this has the same structure as the asymptotic expansions of the one and three 
dimensional systems with short-range potential V, which have no zero-energy eigenstate 
but zero-energy resonance. The existence of the former implies that there is a zero- 
energy eigenfunction belonging to the L 2 -space, while that of the latter corresponds to 
the situation in which there is a function ijj , that is not in I/ 2 (R) but satisfies GqViPq = 
for one dimension or (1 + G V )ipo — for three dimension. The system with the zero- 
energy resonance is known not to be "generic" [ITB, O, 114, I~5|, and in such a case 



(po i n (|3.4j ) is replaced by ipo. In this sense, the one-dimensional free particle system 
is considered to be exceptional. On the contrary, the three-dimensional free particle 
system seems to be generic, since the t _1 / 2 -term does not appear in the expansion of 
exp(— UHq). This is because the free resolvent for the three dimensional case has no 
singularity at the origin, while it appears in ( [2.1| ) (see also [0]). To be precise, the 
asymptotic expansion of the free resolvent for the three-dimensional case is 



Ro(z) = »' 2 K/y 



(3) 



3=0 



(3) 

where Gj is the integral operator with the kernel |x 



yp _1 /47rj! and x, y G R (14 



Therefore the asymptote of exp(- 

r 3/2 



UHq) for the three-dimensional case becomes 

cxv(-m t ) = ••(!/) :i2 (^ 3) ,- )^ 3) +o(r 3/2 ), t 



oo, (3.5) 
3 / 2 is the zero-energy 



in B(s, s'), with large enough s and s'. Here (p ; (x) := (2n 
eigenfunction of H . Notice that does not yield the t _1//2 -term in the expansion 
series, unlikely in the one-dimensional case. 

The formula ( |3.3|) does not bring us the information at each point x. However, it 
is useful for calculating the quantities through the norm or the inner product, such as 
the survival probability. Suppose s, s' > j + 1/2, a 6 R, and M e B(— s',a). Then, 
since Gj is considered as a vector in B(s, — s'), we see that HMGj^Ho- is well defined for 
all ip G L 2 ' S (R). Therefore, for s, s' > max{3n + 3/2, 5/2} and M e B( 
from (U) 



we have 



\M exp(-itH )iJ)\\ v - n 



£(_iy-i r (j + 1/2) (it)-i-Va MG 2 jip 



j=0 



\M\ 



< 



exp(-i^o) - 7T- 1 Y^i-iy-'Tij + l/2)(it)-^ 1 / 2 G 



2/ 



o(t 



-n-l/2\ 



(3.6) 
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for all ip G L 2,S (R). For example, we can take (1 + x 2 ) s '/ 2 E(B) for M with a = —s', 
where E(B) is the spectral measure of the position operator and B an arbitrary 
bounded interval of R, i.e. (E(B)ip)(x) = ip(x) (x G 5) or (i ^ B). Then, 
||Mexp(-i£#o)^ll- s ' = \\E(B)exp(-itH )ilj\\ 2 = J B \ip(x, t)\ 2 dx } and the last quantity 
is called the nonescape probability, which is the probability to find the particle in B at 
a time t. From a similar argument, we also have, 



< 



, exp(-itF )V>> - 7T- 1 E(-l) i_1 r(j + l/2)(itr- 1/2 

3=0 

n 

ex V (-UH ) - 7T- 1 ^(-ly-'TU + l/2)(it)^'- 1 /2 q 2 . 

3=0 



(<f>,G 2j i>) 



o(t 



-n-1/2 



(3.7) 



for s,s' > max{3n + 3/2,5/2} with s > s' and all 0, ip G L 2 ' S (R). The asymptotic 
formula for the survival amplitude of ip is the special case of 



4. Dependence on the initial momentum distribution 

In practical situations, it sometimes happens that there are no contributions from some 
of the G2/S to such quantities like ||Mexp(— itHojipW^ or (0, exp(— itH )ip), when they 
act on a certain vector ip. In this section, we confine ourselves to such situations for the 
survival amplitude {ip, exp(— itH )ip). 

Lemma 4.1 : Let n = 0, 1, If s > n + 1/2 and ip G L 2 ' S (R), then x J ip(x) G 

L X (R) ; for all j = 0,1,..., n. 

Proof : The statement is obtained straightforwardly: f- R \x : 'ip{x)\dx = / R (l + 
x 2 )-( s - j V 2 (l + x 2 y s - j V 2 \xii>(x)\dx < [J R (l + x 2 )-^')dx] 1 / 2 [/ R (l + x 2 )^(x)| 2 dx] 1 / 2 < 
00 for all j = 0, 1, . . . , n with s > n + 1/2. □ 

Lemma 4.2 : Let m — 1, 2, J/s > max{2(m — 1), m} + 1/2 and t/> G L 2,s (R) ; 

t/ien the following three statements are equivalent: 
(a) ip(k) = 0(k m ), k^O. 

(b) ^,G 2 ^)=0, j = 0,l,...,m-l. 

(c) f R x j ip(x)dx = 0, j — 0, 1, . . . , m — 1. 

In particular, we have G^jip =0, j = 0, 1, . . . , (m — l)/2 for odd m, or G^jip = 0, j = 
0, 1, . . . , m/2 — 1 for even m. 

Proof : Suppose that ip G L 2,S (R) with s > max{2(m — l),m} + 1/2. Then by 
Lemma [O] , we first see that x^ip{x) G -^ 1 (R) for j = 0, 1, . . . , max{2(m — 1), m}. Since 
g-ite = ^™-i(_ifc x )i/j! + (_iA ;x ) m / 1 t m - 1 e- ifc:E ( 1 -*)dt/(m-l)!, the fact that (c) implies 
(a) immediately follows, by using 



\m\ 



1 



'2tx Jb. 

< 



\p (— (— ifcr) m ^ ^-^-^(l-i)^ 
.3=0 

h\ m 



jl [m — 1)! jo 

x m ^(x)\dx = 0{k m ), k 



ip{x)dx 



(4.1) 
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To prove the fact that (a) implies (c), let us remember that if)(k) = 0(k m ) (k — > 0) 
means that for some 5 and some finite C > 0, \ip(k)/k m \ < C for all k satisfying 
< \k\ < 5. Then, from Q4.1|) , J R x^if)(x)dx should vanish for all j = 0, 1, . . . , m. The 
fact that (c) implies (b) follows straightforwardly from the identity, 

-2(2j)\^,G 2j ^) = J2(-l) 2 ^( 23 ) ( x 2 ^W)dx f yty(y)dy, (4.2) 



i=o \ " / " R JR 

where the bar (~) denotes the complex conjugate. To prove the fact that (b) implies 
(c), we first use the assumption that —2(ip,G ip) = \ f R ip(x)dx\ 2 = 0. Then from (|4.2j ) 
for j = 1 we have | f- R xip(x)dx\ 2 = 0, from which the remaining equalities recursively 
follow. For the proof of the last part of the lemma, we note that for s,s' > 2j + 1/2 and 

if) e £ 2 ' S (R), 

G 2j ip = <=> \\G 2j il)\\-s' = & f x i ijj{x)dx = (z = 0, l,...,2j). 

Hence, if (c) holds, we obtain the equality G 2 jip = for all j satisfying 2j < m — 1. 
This completes the proof. □ 

Now we shall derive the asymptotic formula for the survival amplitude by combining 



Theorem qJ] with Lemma [O]. The asymptotic formula itself immediately follows 
from (|3.7f) with (j) = if>, under the assumption in Theorem [3.1| . We also see that 
max{3n + 3/2,5/2} > 3n + 3/2 > max{2(m — l),m} + 1/2 for m < n, and thus 



the assumption in Lemma [12] is included in that in Theorem |3.1| . Hence, we finally 
obtain the following theorem for the survival amplitude of if), which is closely connected 
to the behaviour of if) at zero momentum. 

Theorem 4.3 : Let m,n,m < n be non-negative integers, s > max{3n + 3/2, 5/2}, 
and if) G L 2 ' S (R). Then it follows that 

n 

<V,exp(-itff )V>> = n-^i-iy-^U + l/2)(it)^- 1 / 2 (if;, G 2j iP) + o^' 1 ' 2 ). (4.3) 

3=0 

In particular, for some m > 1, if ip(k) = 0(k m ) as k — > 0, then (ip,G 2 jif)) = for all 
j = 0, 1, . . . , m — 1, and vice versa. 



The asymptotic formula (4.3) can also be written as the form without use of G 2 /s. 
Let n = 0,1,.... If s > n + 1/2 and if) G L 2,S (R), ip{k) is n-times continuously 
differentiable in k with if)( n '(k) = (— i) n (27r) -1 / 2 f K x n e~ lkx if)(x)dx. Then, we have from 

(r.a>,„r) = -—— -A - | / x m iP(x)dx = -izIl^M (0 )| 2 . (4.4) 



Lemma |J and (|J) that if ip(k) = 0(k m ) as k -> 

(-l) m (2m 



2{2m)\\m J Jn (ml) 2 

This expression also holds for m = [see ( |3.4| )1. Therefore, we obtain from ( |4.3| ) the 
asymptotic formula for the survival probability 

|(^,exp(-i^ )^)| 2 = r 2 ™- i r(m +Y 2)2 ^(™>(0)| 4 + oft- 2 - 1 ). (4.5) 

(ml r 
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In Theorem |4.3| , the assumption that ip G L 2,S (R) with sufficiently large s is technically 
required. According to the expression in ( |4.5| ), it is worth reviewing this assumption 
in the momentum representation. Let ip G L 2 ' S (R) (s > 0) and [s] denote the smallest 
integer less than or equal to s. Then, by the Plancherel theorem, we have 

oo> f (l + x 2 ) n |^(x)| 2 dx > f \x n 4>(x)\ 2 dx = f \^ n) (k)\ 2 dk 



R, J R, J R 



for all n = 0, 1, . . . [s]. It should be noted here that ij)(k) is implicitly guaranteed to be 
[s]-times differentiable. Hence, as an obvious case, we can find the following subspace 

V := {ij G C°°(R k ) | G L 2 (R k ), n = 0, 1, . . .}, (4.6) 

which satisfies T> C {ip G L 2 (R k ) | if) G L 2 ' S (R)} for all s > 0. Then, Equations (Q and 
( ^7|) with an arbitrary n can be applied to the wave functions belonging to the above 
T>. Examples of such (initial) wave functions include k l /(l + k 2m ) a and k n exp(— k 2 ), 
where I, m, n — 0, 1, . . . and a > with 2ma — I > 1/2. 

In order to see some implications of Theorem |4.3| , let us refer to the following 
two examples. We first consider the rapidly decreasing functions, 4> m {k) = 
N m k m exp(— aok 2 ), as initial wave functions, where m = 0, 1, . . ., ao > 0, and N m := 
[r(m + l/2)/(2ao) m+1 ^ 2 ] -1 ^ 2 being the normalization constants. Then it is obtained 
through the Laplace transform 

/oo 
\4> m (k)\ 2 exp(-itk 2 )dk 

= [1 + it/(2a )]-( m+1 / 2 ) = (it/2a )- (m+1/2) [1 + Oft" 1 )]. (4.7) 

On the other hand, we see that the right-hand side of ( [4.5| ) exactly corresponds to that 
in ( |4.7| ) in the leading order. The other example is a special case that V»(0) = for 
all n = 0, 1, . . .. It is worth noticing that for such a initial wave function we clearly see 
from (4J5) that 

|(^,exp(-itW)r = o(t- 2n - 1 ), 

for every n > 0. That is, the survival probability decays faster than any power of 
t -1 . However, it must decay slower than any exponential at long times for Ho > 
]l||. This strange decay behaviour is also found in a study of the time operator 
(Proposition 3.2 in |TJ]). The set of such a special wave function is given, e.g., by 
d := {4) G C °°(R fc ) | 3k > 0;i)(k) = 0, for k G [-k ,k Q }}. We see from (p|) that 
C\ C T>. A wave function in C ; has a positive lower-bound /c 2 on energy. For instance, 
the following function 

2 (h) _ / exp(-l/[A; 2 - (k - d) 2 }) (\k - d\ < k ) 
m ~\ (\k-d\>k ) ' 

where d > k > 0, surely belongs to C[. 
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5. Concluding remarks 

We have derived the asymptotic expansion of the time evolution operator for the 
one-dimensional free particle system, in terms of the operators which are expansion 
coefficients of the free resolvent at small energies. This enables us to obtain the 
asymptotic formula for the survival probability of ip, and also to evaluate, in a systematic 
way, the condition for the initial wave function if> which makes the first several terms 
of the asymptotic formula vanish. We have found that if ip(k) = 0(k m ) for some 
non-negative integer m at zero momentum, the asymptotic power of £ _1 for the survival 
probability must be 2m+l. In other words, the information about the initial momentum 
distribution in the vicinity of zero momentum is reflected in the asymptotic decay form 
£-2m-i a £ i ori g times. Our results are essentially due to the choice of the initial wave 
functions ip in L 2 ' S (R) with sufficiently large s. This guarantees the existence of the 
higher derivatives at zero momentum (see Remark ??). However, there is another wave 
function such that ^^(O) = up to n = m — 1, while its m-th derivative ^ m \ti) 
diverges. Related wave functions are considered in || ||, [L(| [Tl| . For such states, the 



asymptotic formula ( |4.5|) is not correct. Indeed, the actual asymptotic decay form of 
the survival probability includes terms of non-odd power of t _1 . We hope to address 
this issue in the future. 

Acknowledgements 

The author would like to thank Professor I. Ohba and Professor H. Nakazato for useful 
and helpful discussions. 

Appendix A 

Dependence of the survival probability on the initial state ip f° r one- dimensional free 
particle system can be seen in the following inequality (Theorem 4.1 and section VI in 

0), 

4||T ^ 2 |H| 2 
t 2 



|(^exp(-itifoMr< " °Z > teR ' ( A1 ) 



where To is the Aharonov-Bohm time operator [Tj| . This brings us with an interpretation 
of ||T ^||, or of the time uncertainty calculated from T . The relevant information 
herein is that for any L 2 -function tp whose several moments are finite, e.g., ip G L 2,S (R) 
(s > 3/2), we have 

||T VH| < oo <=> ijj(k) = 0(k 2 ), k -> 0. (A2) 



This implies, together with (|A1|) , that the condition (|A2|) at zero momentum imposes 
on the survival amplitude (tp, exp(— itiJo) 1 ^) a decay faster than t , which is obviously 
faster than t" 1 / 2 of the usual decay law for the one dimensional free particle system. 
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To prove the relation ( |A2| ) , we first define the Aharonov-Bohm time operator T , 
which is mathematically well treated in the scheme of the axiomatic quantum mechanics 
fl2|, |20|. We define this operator as follows: the domain of T is 



D(T ) := |^ G L 2 (R) lim ^(k)/\k\ 1/2 = and (T i/;)(k) G L 2 (R k ) } , (A3) 



and its action 



i 



fdi>(k)/k ld^(fc)\ 



W)(k) = - I T dk > + ^-^J , a.e. k G Rfc, ip G D(T ), (A4) 

where ip(k) is assumed to be different iable everywhere except the origin. For a ip that 
belongs to L 2,S (R) (s > 3/2), ip(x) and xi[)(x) in -^ 1 (R) by Lemma O , This implies 



^ip(k) to be differentiable everywhere including the origin. Then, it follows as in (4.1) 
that for k — > 

V>(fc) = ^(0) + fc*0'(O) + 0(k 2 ), $'(k)=$'(0) + O(k). (A5) 

Hence 

for small fc. Furthermore, since ^ G L 2 ' S (R) C L 2,1 (R), J- R \xip(x)\ 2 dx = 
J Rk \ip'(k)\ 2 dk < oo. Thus, (T ^)(k) is assured to be square integrable on (— oo, —6} U 
[5, oo) for an arbitrary 5 > 0. 

Now, to prove the fact that ||T ^|| < oo implies ip(k) = 0(k 2 ), k — > in 
the relation flA"2j), let us suppose that (To^)(fe) G L 2 (Hk)- Then, from the general 
property of L 2 -functions, we see that (Toip)(k) G L 1 ([— 5, 5]). This contradicts with 
( |A6| ) unless tp(0) = ^'(0) = 0. Thus we have that ip(k) = 0(k 2 ). Conversely, if 
i)(k) = 0(k 2 ), it follows from flX5]) that ^(0) = $'(0) = 0. Then, (|A|) implies that 
(T xp)(k) G L 2 {[— 5, 6]). Hence, (T ■?/>)(&;) belongs to L 2 (R fc ), and the proof of the relation 
is completed. 



Appendix B 



In this appendix, we derive the formula ( |3.2j ) which directly relates exp(— UHq) to the 
spectral density. Let us remember that the time evolution operator for one dimensional 
free particle system is explicitly represented as 

ip(x,t) = (exp{-itH )il>)(x) = (4vrifr 1/2 f exp(i|x - y\ 2 /At)i)(y) dy = 0(r l/2 ) 



R 



(Bl) 



for all i) G L 1 (R)HL 2 (R) and t ^ 0. Note that L 2 ' S (K) C L*(R) for s > 1/2. Then, 
exp(— itH ) is considered as an integral operator belonging to B(s, —s') for s, s' > 1/2. 
Let s, s' > 1/2 and ip G L 2 ' S (R). We have an equality 



r-R 



E'(X)e-' ltx dX^ 



R 



dy 



R cos(X 1 ^ 2 \x — y\)e~ 



ifA 



27TAV2 



-dA 



ip(y), a.e. x G R, 



The power decay of the survival probability at long times 
for positive r and R. Then it follows from ( |B1| ) that 

J* E'(\)e- itx d\ilj - exp(-itH )iP 

« cos(X^ 2 \x - y\)e~ itx 
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< 



R JR 



27TAV2 



dA 



exp(— \x — y| 2 /4it) 



(47ri*)V2 



dxdy 



(l+x 2 ) s '(l + y 2 y' 



We here see that 



sup 

x£TL,R>0 



r 



cos(A 1 / 2 |x|)e- itA 



-dA 



< oo. 



(B2) 



(B3) 



10 2ttA 1 /2 

To derive this, we have assumed t > 0, however, the following argument can be applied 
to negative t, 



^cos(A 1 / 2 |a;|)e- 

27TAV2 



if A 



-dA 



(27T) 



-1 



(2 



7T 



v-1 



|^ /2 exp(-it(e + a) 2 )de 
/ eX p(-i^ 2 )de+ / /9 exp(-it£ 2 )d£ 

JO J-R 1 / 2 +a 



where a = \x\/2t. The last two integrals are estimated, by the complex integrals with 
the contour in fourth quadrant, to be 



exp(-it£ 2 )d£ 



Trfl-expf-tl^ial 2 )) r\R x ' 2 M 

n 1 — exp(— tx 2 ) 1 fW 

< — sup 1 — \ — < oo. 

~ 4t x>o x 2 V t 



This leads to (p3|). We also see straightforwardly 

f R : cos(A 1/2 |x - y\)e~ itx „ exp(-|x - w| 2 /4it) 
hm hm / — dA = — r-rjz . 

Therefore, by the dominated convergence theorem, we have from (|B2|) , ( p3|) , and 



(B4) 



sup 

V>eL 2 > s (R),V>^o 



exp(-itH )i/i - f E'(X)e-' ltx d\ij 

Jr 

as r j. and R —>■ oo, and we finally obtain the formula (|3.2j ). 



Appendix C 



In order to prove Theorem |3.1| , we first summarize the several property of E'(\). By 
Lemma |2.3| , we see that for n > the remainder F n (X) in is (n + l)-times 

differentiable, in B(s, — s') for s,s' > 2n + (n + 1) + 1/2 = 3n + 3/2, and satisfies 
that (d/dA) r F„(A) = o(A"- r - 1/2 ) as A | (r = 0, 1, . . . , n + 1). On the other hand, 
we see from Lemma that (d/dA) r £'(A) = 0(A^ r+1 )/ 2 ) as A — > oo, in B(s, -s') 
for s,s' > (m + 1) + 1/2 = m + 3/2 (r = 0, 1, . . . , m + 1). In particular, if m > 1, 
(d/dA) m+1 .E'(A) is integrable on [S, oo) for an arbitrary 5 > 0. 

Let us now split the integral in ( |3.2|) into two parts by writing 

e'(A) = 0(A) j e;'(A) + (i-0(A)) j e;'(A), 
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where G Cq°([0, oo)) and satisfies 0(A) = 1 in a neighbourhood of A = 0. Such 
a function is realized by /(A) = 1 — J A g(x)dx, where g(x) = h(x)/ J- R h(x)dx and 
h(x) = exp(-l/[l - (re - cf) 2 ]) (|x - d| < 1) or (|x - d| > 1) with d > 1. 

From Lemma 10.1 in Jl4| and the discussion as mentioned above, we see that 
(1 — 0(A))£"(A) has a contributions of o(t~ m ~ l ) to exp(— itffo) in B(s, —s'), where m > 1 
and s, s' > m + 3/2. 

On the other hand, the contribution of 4>(X)E'(\) to exp(— itH ) gives the main 
part of the asymptotic expansion. Then, the coefficient of Gij is given by 

I" <P(X)X'- l ' 2 e- itx dX =i j ^-\ r A~ 1/2 e- ia dA + /°°(0(A) - l)A~ 1/2 e~ iu dA 
Jo dP Uo Jo 

= r(j + l/2)(it)-^ 1 / 2 + ^ — j o (0(A) - l)A-% itt dA. (CI) 

Note that since 0(A) — 1 and all its derivatives vanish in the neighbourhood of A = 
and G C£°([0, oo)), the last term in ( |Cl| ) decays faster than any negative-power 
of t. Furthermore, we understand, from Lemma 10.2 in and the discussion in 
the first of the appendix, that if s,s' > 3n + 3/2, any contribution of the Fourier 
transform of the remainder 0(A)F n (A) to exp 

(-it-ffo) is of oit- 71 - 1 / 2 ) in the norm of 
B(s, —s'). Summarizing the above arguments, we finally obtain, under the condition 
s, s' > max{3n + 3/2,5/2}, 



< 



exp(— itH ) — 
exp(-itHo) 



71 



3=0 



E\X)e~ ltx dX 



+ 



/ E'(X)e~ itx dX - EC-l^rO' + l/2)(it)^'- 1/2 G 2j 

Jr 3=0 



+ o(t 



-n-l/2\ 



as r I and R — >• oo. This is just the asymptotic expansion of exp(— itH ) in (|3 
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